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Abstract. A new method is developed to calculate the amplitude of the galactic warps generated by a torque due 
to external forces. This takes into account that the warp is produced as a reorientation of the different rings which 
constitute the disc in order to compensate the differential precession generated by the external force, yielding a 
uniform asymptotic precession for all rings. 

Application of this method to gravitational tidal forces in the Milky Way due to the Magellanic Clouds leads to a 
very low amplitude of the warp, as has been inferred in previous studies; so, tidal forces are unlikely to generate 
warps, at least in the Milky Way. If the force were due to an extragalactic magnetic field, its intensity would have 
to be very high, greater than 1 fj,G, to generate the observed warps. 

An alternative hypothesis is explored: the accretion of the intergalactic medium over the disk. A cup-shaped dis- 
tortion is expected, due to the transmission of the linear momentum; but, this effect is small and the predominant 
effect turns out to be the transmission of angular momentum, i.e. a torque giving an integral-sign shape warp. The 
torque produced by a flow of velocity ~ 100 km/s and baryon density ~ 10 -25 kg/m 3 is enough to generate the 
observed warps and this mechanism offers quite a plausible explanation. First, because this order of accretion rate 
is inferred from other processes observed in the Galaxy, notably its chemical evolution. The inferred rate of infall 
of matter, ~ 1 Mo/yr, to the Galactic disc that this theory predicts agrees with the quantitative predictions of 
this chemical evolution resolving key issues, notably the G-dwarf problem. Second, because the required density 
of the intergalactic medium is within the range of values compatible with observation. By this mechanism, we can 
explain the warp phenomenon in terms of intergalactic accretion flows onto the disk of the galaxy. 

Key words, intergalactic medium — galaxies: structure — Galaxy: structure — galaxies: interactions — galaxies: 
kinematic and dynamics — galaxies: magnetic fields 



1. Introduction 

Many spiral galaxies present warps, distortions to a flat 
disc with an integral-sign shape. The Milky Way is an 
example (Burton 1988; 1992). Indeed, most of the spi- 
ral galaxies for which we have relevant information on 
their structure (because they are edge on and they are 
nearby) present warps. Sanchez-Saavedra, Battaner & 
Florido (1990) and Reshetnikov & Combes (1998) show 
that nearly half of the spiral galaxies of selected samples 
are warped, and many of the rest might also be warped 
since warps in galaxies with low inclination are difficult to 
detect. 

The intergalactic magnetic field has been suggested as 
the cause of galactic warps (Battaner et al. 1990; Battaner 
et al. 1991; Battaner & Jimenez- Vicente 1998). This is in 
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our opinion a serious proposal (an opinion not held, how- 
ever, by Binney 2000) which could explain many of the 
observations, although observational support is still con- 
troversial. The postulated alignment of warps of different 
galaxies (Battaner et al. 1991) and the differences between 
the gaseous and stellar warps (Porcel et al. 1997) can have 
alternative explanations as we shall see in the present pa- 
per. 

Gravitational tidal effects on the Milky Way due to the 
Magellanic Clouds are not enough to justify the observed 
amplitude of the warp. Hunter & Toomre (1969) calcu- 
lated that the Clouds with mass M sat = 10 10 M and dis- 
tance d — 55 kpc would generate a warp of amplitude less 
than 117 pc at radius 16 kpc in the most favourable case 
instead of the observed 2 or 3 kpc. The Magellanic Clouds 
are near the pericentres of their orbits around Milky Way 
(Murai & Fujimoto 1980; Lin & Lynden-Bell 1982), so it 
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is not expected that this amplitude could be greater due 
to a closer approach of the Magellanic Clouds in a recent 
past. 

Weinberg (1998) proposed a mechanism to amplify the 
tidal effects due to a satellite by means of an intermediate 
massive halo around the galactic disc, but Garcfa-Ruiz et 
al. (2000) have found that the orientation of the warp is 
not compatible with the generation of warps by means of 
this mechanism if the satellites are the Magellanic Clouds. 
Quantitatively, a better prospect would be the Sagittarius 
dwarf galaxy (Ibata & Razoumov 1998) since tidal effects 
are proportional to ^ at , the galactocentric distance to 

"■sat 

this dwarf Galaxy is only 16 kpc and its mass ~ 10 9 M Q 
(see Q. 

Binncy's (1992) review concludes that halos dominate 
the dynamics of the warps, although he also points out 
that "warps will in the end prove to be valuable probes of 
cosmic infall and galaxy formation". In a subsequent pa- 
per (Jiang & Binney 1999), cosmic infall is used to explain 
the reorientation of a massive Galactic halo (8 degrees per 
Gyr) which produces a warp in the disc. This model re- 
quires a halo ten times more massive than the disc, an 
extremely high accretion rate (3 disc masses in 0.9 Gyr) 
and, in this scenario, after a sufficiently long time, the 
angular momentum of the Galaxy would become parallel 
to the direction of the falling matter causing the warp to 
decay. This last problem might be solved by including a 
prolate halo (Ideta et al. 2000). The general case of warps 
produced by the dynamical friction between a misaligned 
rotating halo and disk was also studied by Debattista & 
Sellwood (1999). Other proposals which include a massive 
halo have also serious flaws or not very plausible assump- 
tions (Nelson & Tremaine 1995; Binney et al. 1998). 

In spite of the importance which is given to the halo 
in the dynamics of galaxies, it is quite possible that they 
may not even play a major role in the formation of warps. 
The mass of the halo of the Galaxy is not especially well 
determined and the mass fraction of galaxy halos might be 
small (see, for instance: Nelson 1988; Battaner et al. 1992; 
Evans 2001 and references therein). We should emphasize 
here that the presence or absence of a very massive halo 
will not modify qualitatively the arguments presented be- 
low and would imply quantitative changes within the same 
order of magnitude. We will argue in this paper that warps 
can be generated without massive halos or magnetic fields, 
although our results are perfectly compatible with the ex- 
istence of these. 

We aim here an alternative solution to those hith- 
erto hypothesized. The mechanism of generation of warps 
which is explained in this paper solves these previous diffi- 
culties and does not require implausible assumptions. It is 
even simpler than the hypotheses previously proposed; it 
requires only the infall of a very low density intergalactic 
medium onto the disc without the dynamical intervention 
of an intermediate halo. It is a very simple idea but it 
works well, as will be shown below. 



An analytical calculation is performed to reduce the 
problem to a differential equation and some integrals, 
which are subsequently solved by means of numerical algo- 
rithms. A new method is developed in §|| to calculate the 
warp parameters from an induced external torque which 
takes into account the interaction between all rings of the 
galactic disc. The external torque induced by the accre- 
tion is calculated in §||, which allows us to estimate the 
required density of the inflow, in §|j. 

2. Warp induced by an external torque 

In this section, we will explain the mechanism of gener- 
ation of warps in a galactic disc due to a net external 
torque. This general method is applicable to any kind of 
torque acting over the disc and will be used to derive the 
properties of the warp induced by an intergalactic accre- 
tion flow. 

In order to describe a warp in a galactic disc, we will 
use the usual model of tilted rings (Rogstad et al. 1974): 
the disc is taken to be a set of concentric rings, each of 
radius between Ri and Ri + dR having angle inclination 
on with respect to the central disc and intersecting its 
plane at two nodes, perpendicular to the points on the 
ring where the elevation is maximum and minimum. The 
line joining the nodes ("line of nodes"), will be taken as 
common to all the rings, i.e., the nodes for all rings are 
aligned, as observed in our Galaxy; this is the equilibrium 
state in our model. Therefore, the only parameters which 
define the warp are: the direction of the line of nodes and 
a function, a(R), the maximum angular elevation of the 
ring of radius R with respect to the plane defined by the 
central disc. 

A torque applied to a rotating rigid body (in our case, 
the rigid body is a ring) produces a precession in it, as 
in the case of the equinox precession of the Earth. The 
importance of precession in galaxies was indeed first rec- 
ognized by Lynden-Bell (1965). This section describes how 
the differential precession in successive rings generates a 
warp due to their redistribution. The general formalism 
included here is applicable to any kind of torque induced 
by a external force over a set of nearly coplanar rings. 
In the following section, the calculation of the torque for 
the specific case of the intergalactic accretion flows will be 
derived. 

Other analytical approaches to the problem were used 
previously. For instance, the consideration of the warp as a 
product of equilibrium forces of vertical components, such 
that F oxt , z = -Fg raViZ , where F grav is the gravitational force 
due to an axisymmetric potential (used in Kahn & Woltjer 
1959; Binney 1991; Binney 1992 (his section 2: "Naive the- 
ory"); Battaner & Jimenez- Vicente 1998). Since the orbit 
within a ring is tilted, the centrifugal forces will cancel 
both radial and vertical components of the gravitational 
forces so this consideration is not valid. Our approach has 
some advantages and finer calculations than these papers. 
A mass in the centre of the galaxy cannot provide any 
torque to the rings. However, the axisymmetric poten- 
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tial of the disc is distorted by the warp (Binney 1992), 
and these non-axisymmctries are those responsible for the 
warp itself. The differential orientation of the successive 
rings rather than a point mass placed in the centre pro- 
duces a torque. The approach used in Hunter & Toomre 
(1969) is much better although the analysis is very differ- 
ent to the one presented here. 

2.1. Equations 

We consider the disc to be made of material on circular or- 
bits with angular velocity w ro t (R) ■ If initially all the orbits 
are in the same plane (the plane of the disc), perpendicu- 
lar to a unit vector k, under the influence of an external 
torque with a non-vanishing component perpendicular to 
k the orbits will precess. If the precession is not equal 
for all orbits they will not remain in the same plane, so 
that the vector perpendicular to the plane of a given orbit 
k(R, t) will be a function of R and of time t. 

In the general case, the component of the external 
torque along the vector k(R,t) will induce changes in 
w r ot(-R)- I n addition, other moments of the external force 
may produce changes in the shape of the orbits. However 
in the applications relevant to the present problem these 
changes are negligible and we will take the orbits to remain 
circular, and with essentially constant Lo mt (R). In this case 
the dynamics of the disc under an external torque may be 
reduced to an equation in partial derivatives for k(R,t), 
or in other words a set of an infinite number of ordinary 
differential equations, (one for each value of R). The unit 
vector k(R,t) is defined by two parameters which are in 
fact angles. 

The equations for the evolving system, (throughout 
this work we use equations valid for an inertial, non- 
rotating frame, and at no stage do we use a rotating frame) 
averaged over times longer than the orbital periods are: 



dJ_ 

dt 



(R,t) =T(R,t) = To* (ii,t)+ T int (#,«), 



(1) 



where J(R, t) is the angular momentum for the ring with 
radius between R and R + dR, and T ox t, T; n t represent the 
average torque in the ring due to the external force, and 
to the gravitational interaction with the rest of the disc 
respectively. 



2.1.1. External torque 

The external torque is: 

Text(-R^) = T ext [k(R,t).u,R]i 



(2) 



io = jo x k , 



where it is a unit vector in the direction associated with 
the cause of the external torque, which may be e.g., the 
direction of the infalling gas flow or the position of a per- 
turbing galaxy, j is a unit vector parallel to the pro- 
jection of u onto the initial plane (perpendicular to fc ), 



T oxt (x, R) is a function of the cosine of the angle between 
it and k(R, t), which characterizes the specific mechanism 
giving rise to the torque. 

2.1.2. Internal torque 

The internal torque is: 



Ti„tOR,t) = / [dr s (R,t)dR]dS, 

J all rings 



(3) 



and drs(R,t)dR dS is the gravitational torque that the 
ring of radius between S and S + dS produces on the ring 
of radius between R and R + dR, which is: 



dr s (R,t)dR dS = Ga(R)a(S)dR dS 
XT int (k(R,t).k(S,t), R, S)k(R,t) x k(S,t), 



(4) 



where Tint (cos ctR.s, R, S) is a well determined function 
(the same for any mechanism) of the cosine of the angle 
a_R,s between the planes of the two orbits. Explicitly: 



= Ga(R)a(S)dR dS 



q2 rZ-n [2 



d(f>2 



(5) 



[1 + {S/R) 2 - 2{S/R){sin4>i sin0 2 + cos0i cos0 2 cosa fl , s )] 3 / 2 
x [sin 0i sin «ij,s cos 4>2t{R, S, t) — cos 0i sin <xr } s cos 4>2j{R, S, t) 



+ [cos 4>i sin 02 — cos 02 sin0i cos oiR z s]k(R, t)], 
a R>s = a{S, t) - a{R, t); a(0, t) = 0; 
cosa.R,s = k(R,t).k(S, i) 



i(R, S, t) ee j(R, S, t) x k(R, t); j(R, s, t) = 



k(R,t) x k(S,t) 
\k(R,t) x k(S,t)\ 



The torque between two rings, drs(R, t)dR dS, is pro- 
portional to —j(R, S, t), and the external torque must be 
parallel to this. The z-component of the external torque, 
if this were present, would produce an acceleration in the 
rotation of the disc rather than warping. 

Other components of a galaxy produce negligible con- 
tribution. The bulge in practice contribute negligibly to 
the torque: firstly, because it is more spherical than the 
rings and a spherical distribution of mass produces no 
torque; and, secondly, because the distance of the bulge 
to the outer rings is large enough for these to produce 
negligible effects (the torque is proportional to S 2 /R for 
small S, where S is the radius of the structure). Numerical 
experiments were carried out which confirmed this point. 
A massive halo, if it exists, would produce an extra in- 
ternal torque if the halo is non spherical. In practice, a 
massive non-spherical halo will change quantitatively the 
amplitude of the warp but, qualitatively the mechanism 
will be the same; the warp amplitude will be reduced, 
never increased, because a massive halo would keep the 
rings more tightly bound. At present, no precise calcula- 
tions including non-sphcrical massive halos are given here, 
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but the net result would be to increase the required ac- 
cretion rates inferred below. Only the inner mass of the 
halo ellipsoid with semimajor axis equal to R produce net 
torque, assuming there is a constant ellipticity halo, whose 



quadrupole is 



4 vLttR fh(R) 



{v ro t (R) is the rotation ve- 



15 ~G T " 

locity of the Galaxy; e is the eccentricity of the ellipsoids; a 
mass distribution for the halo derived from a hypothetical 
flat rotation curves due to a dark halo is assumed; fh(R) 
fraction of the mass M(R) embedded in the halo) while 
the disc quadrupole component is ~ 12T7cr(RQ)e Ro ^ hR hj i 
for a model such as (jl^). The disc dominates at R < 22 
kpc for e = 0.2 and f h = 0.5 (Kuijken & Dubinski 1995), 
and the contribution of the halo for R < 16 Kpc is less 
than 40% of the disc contribution. This means that the or- 
der of magnitude of the torque will be not affected by the 
inclusion of the massive halo, and will be perhaps affected 
by a factor less than ~ 1.4 at R < 16 kpc. 

For small angles, in a linear approximation, the pro- 
portionality which results is: 



lim d,Ts{R)dR dS ct —an.sj{R, S,t). 

ob.s >0 



2.1.3. Dynamics and evolution of the warp 



(G) 



The precession velocity of k(R,t) is much less than 
Wrot t) , so we have 



J(R,t) m Iuj Iot (R,t)k(R,t) 
I = 2irR 3 a(R)dR, 



(7) 

(8) 



where I is the moment of inertia of the ring, and er(i?) the 
surface density of the disc. 



a(R) 



dzp disc (R, z), 



(9) 



where pdisc is its spatial density, independent of <f> in an 
assumed axisymmetric case. 

This approximation is typical in planetary precession, 
for instance. Using this expression in (]l|), and calling ry, 
t ±_ the components of r along k(R,t) and perpendicular 
to k(R,t) respectively, we obtain 

\Icb Iot (R,t)\ = |T||| = \k(R,t).r(R,t)\; (10) 

dk(R,t) = r± 
dt ' Iu> rot (R,t) 

_ [r(R,t) - {k(R,t).T(R,t))k(R,t}} 
Iu TO t(R,t) 

With the set of initial conditions k(R,ta) = fco = 
(0,0, 1) Vi?, equations ( p^| ) can be integrated to give the 
configuration of the orbits, and hence the shape of the 
disc, at any time. 

The method described here is not restricted to linear 
perturbations but is valid as long as the deformed disc 
is describable by a set of tilted rings. However, for large 



deformations the orbit cannot remain flat and the approx- 
imation will fail. This is not the case however in the exam- 
ples given in the present discussion. It must be noted that 
equation ( |i"(i| ) for (4jr) is a first order differential equation, 
which implied that it is the angular velocities and not their 
time derivative that are determined by the torque. So no 
matter how large (modulus of vector 4|) may be, if the 
torque vanishes instantly so will 41. This is due to the 
fact that the energy associated with 4Jr is much smaller 
than the energy corresponding to Lu ro t] the full equations 
are obviously second order. 



2.1.4. Equilibrium configuration 

Here we are interested in the stationary, or equilibrium 
configuration that we assume exists, rather than in the 
evolution. It is clear that in this situation all orbits (or 
rings) must precess around u at the same constant speed, 
u>p, keeping k(R,t).u constant in time, though different 
for different orbits. For a(R) and the precession velocity, 
u)p(R), to be independent of time, r(R,t) must be per- 
pendicular to k(R,t) and u for any R. Since this holds 
for T ex t(R,t) it must also hold for T- m t(R,t). However, 
examining the expression for this last vector in (|i"o|), for 
this condition to hold for all orbits, they must all intersect 
along the same straight line, parallel to the unit vector i. 
In this case the position of the orbits is simply given by 
the function a(R) = cos _1 (fc(i?, t).feo). 

If the vector k(R,t) precesses around u with angular 
velocity uj p (R) we then have: 



tOp{R)u x k{R, t) 



dk 
~dt' 



(11) 



From this equation we obtain w p (R)[a(R)] as a func- 
tion of a(R). Now, since one of the conditions for a station- 
ary configuration is that oj p (R) be independent of R, set- 
ting the derivative of this function with respect to R equal 
to zero provides us with the functional equation which we 
must solve for a(R) to give the shape of the distorted disc, 



du>p(R) 
dR 



a(R)] = 0. 



(12) 



The functional equation ( |12|) can be solved numerically 
using the method explained in appendix From this, 
we can obtain a(R), i.e. the amplitude of the warp as a 
function of the distance from the Galactic centre for a 
given external torque r ex t- 

2.1.5. Few remarks on the transient regime towards 
the equilibrium configuration 

As explained above, we are not concerned here about the 
transient regime. However, it is interesting to comment on 
the qualitative aspects of this regime, since they indicate 
that some of the mechanisms considered here may well be 
of interest. For a mechanism to be acceptable it should 
lead to a stationary result, in reasonable agreement with 
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Fig. 1. Spherical triangle which relates different parame- 
ter in differential precession of two rings. 

the observations, which are assumed to correspond to sta- 
tionary systems. However the mechanism is not plausible 
unless it offers a credible scenario for the formation of a 
warp once the external influence initiates its operation. 
We will now see how qualitative considerations relative to 
the warp formation process do indicate different degrees 
of plausibility for different mechanisms, even where their 
predictions for the stationary state are similar. 

To carry through this qualitative analysis, it is useful 
to consider eqns. (10) for a case in which all the mass is 
concentrated in two rings (orbits), one of these being much 
more massive than the other. r ox t is always perpendicular 
to u, k (although this does not hold for Tint), and since 
initially T; nt = the initial motion is a pure precession of 
both rings around u. However this precession occurs, in 
general, at a different rate for each ring. So the plane of 
one of the rings will be elevated "above" the plane of the 
other, at a dihedric angle a. After they have precessed 
differentially through an angle (3, and the line of nodes 
forms an angle SI with the projection of the vector u onto 
the plane of the more massive ring, the relation between 
these three quantities and 9q (where sin#o = fc.it) is given 
implicitly by the spherical triangle in Fig. [I]. Clearly as 
P goes to zero, f2 also goes to zero. This means that the 
initial node direction is at tt/2 with respect to the final 
stationary position. Now as a increases r increases, pro- 
ducing an additional precession of the plane of the lighter 
ring around the vector perpendicular to the plane of the 
heavier. 

If the light ring precesses initially around u faster than 
the heavy ring (this holds in the general case for rings 
with an initial precession velocity faster than the mean), 
the additional precession induced in the heavy ring is in 
the sense opposed to f2. So this additional precession must 
sweep out an angle somewhat greater than tt/2 before it 
reaches the stationary state. If the light ring precesses ini- 



tially more slowly the corresponding angle would be some- 
what less than tt/2. However the implication is that the 
time required to reach the stationary condition is of order 
a quarter of the initial precession period. 

Translating this to the real case, it implies that the 
time for the formation of the warp must be of the order 
of one quarter of the final global precession period. This 
time could be quite long, and a formation procedure which 
leads to a shorter time scale would be more plausible. This 
is the case, for example, where r cxt is absorbed by the gas. 
In this case to simulate the evolution we must consider 
equation (|h]) for two sets of rings, gas rings, which are 
affected by both T cxt and Tj nt , and stellar rings, which are 
affected only by Xj nt . Qualitatively the gas rings will move 
initially as indicated above, giving rise to a Tj nt which 
is perpendicular to r cx t- The stellar rings which initially 
remain in the same plane start to move under the influence 
of Ti n t, with k(R,t) parallel to jo so that the resulting 
nodes are already in their final position. For very large 
radii i?, the rings would have precesed by a small amount 
before the T;„t because the torque of the inner stellar rings, 
which already are warped, dominate that due to the gas. 
Thus, the line of the nodes of these orbits will form some 
angle (going asymptoticaly to 7r/2) with that of the inner 
orbits. The nodes of the outer orbits trace a leading spiral. 
This seems to agree qualitatively with some observations 
(Briggs 1990). 

The gas rings then move rapidly under the much 
stronger Tj nt generated by the stellar rings, and align gen- 
erally with them, although the final warp will be slightly 
different for the gas rings and the stellar rings, so that the 
T ex t experienced by the gas is transferred gravitationally 
almost entirely (but not quite) to the stellar ring at the 
same radius, although probably the relative displacement 
between the two is smaller than the thickness of the disc, 
and will not be considered in the present paper. 

As we pointed out earlier, equation ( |ffj| ) is first order, 
so if when the orbits arrive at their stationary position the 
torque takes their stationary values (which happens when 
all orbits reach their stationary positions simultaneously) 
the line of the nodes of every orbit will stop precessing 
in the mean plane defined by all other orbits, so station- 
arity is achieved. However this do not seems likely and a 
more probable outcome is that the line of the nodes os- 
cillates around the equilibrium position. Perhaps, friction 
could damp the amplitude of these oscillations. There is 
no problem in this oscillations being eliminated without 
dissipation (since they contain no energy) but it seems 
unlikely because of the degree of conspiracy between the 
orbits required. Thus our stationary solution will corre- 
spond to the statistical mean, but at any time there will 
be present minor wiggles (in a(R)) superimposed on it. 

The comment on warp formation presented here is a 
qualitative anticipation of future detailed work and does 
not affect any of the conclusions of the present work. 
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Fig. 2. Graphical representation of the Galactic disc (with 
coordinates xyz) and a point mass at P (with coordinates 
x'y'z'). 

2.2. Example of application: gravitational torque 

Once we have a general formalism in which a warp is gen- 
erated when a torque is produced in a galactic disc, we 
can analyze possible mechanisms to produce that torque. 
A conventional example is the torque generated by gravi- 
tational tidal effects. 

In Fig. ||, we show the Galactic disc centred at O and 
a point mass at P such that OP = dep. In the case of 
the Milky Way disc, the Sun would be situated in the disc 
with R = Rq, <fi = 180°. The unit vector can be expressed 
as 



ep — cos <pp cos ft pi + sm <pp cos ft pj + sin ft p 



d P k. 



(13) 



It is well-known that the torque produced by a point 
mass of mass mp on any axisymmetric body is, consider- 
ing only the quadrupolar term, 



3G m P {h-h) . 

Tgrav S3 — COS ftp Sin ftp 

x(sin0pi — cos <j)pj), 



(14) 



where I\ and I3 are the inertia tensor components (in this 
case for the disc). 

The solution of (|l^) with r ext = r grav by means of the 
numerical method explained in appendix [X] gives a(R) 
for the galactic warp. It will depend on the adopted 
Galactic model. For the Milky Way, a suitable model can 
be adopted as follows: 



The surface density of the disc, (0), is taken as 



a(R) = l a ( R ^ e ~ E5 ~> R <3 R o 
I 0, R>3R 



(15) 



where the local surface density is a(Ro) — 48 M0PC -2 
(Kuijken & Gilmore 1989), the distance to the Galactic 
centre is Ro = 7.9 kpc (Lopez-Corredoira et al. 2000); and 
the scale length hp — 3.5 kpc (Bahcall & Soneira 1980). 
We truncate the exponential disc at 3Rq. The Galactic 
disc undoubtedly extends to larger radii but the effects 
of those outer rings can be considered negligible in the 
present dynamical context. 

The rotation velocity is taken as 



Vrot{R) = Ru ro t(R) 



200 km/s, R < 15 kpc 

200^/15 kpc/R km/s, R > 15 kpc 



(16) 



(Honma & Sofue 1996). Outside the stellar truncation ra- 
dius (~ 15 kpc), a Keplerian law (w rot oc RT' 1 / 2 ) is fol- 
lowed, which implies that we have neglected dark matter 
contributions from larger radii (according to Honma & 
Sofue 1996, there is no essentialy dark matter beyond 15 
kpc). 

For the Large Magellanic Cloud, we adopt ftp = —33°, 
mp — 10 10 M Q and d — 55 kpc, the same values that 
Hunter & Toomre (1969) adopted. 

The precession we obtain is ui p — 2.7 x 10~ 19 rad/s 
and the amplitude of the warp is shown in Figure |3| where 
a(R) w \z\/R. The functional shape of a{R) resembles 
that of observational data, but the amplitude is very dif- 
ferent: a factor of 20 or 30 separates the two curves, in 
agreement with Hunter & Toomre (1969). This means that 
we would require the mass of the Magellanic Clouds to 
be around 2 or 3 x 10 11 M Q to justify the Milky Way 
warp as a product of the gravitational interaction with 
the Clouds. The calculation of the amplitude with the 
Saggitarius dwarf galaxy (mp = 10 9 M , d — 16 kpc; 
Ibata & Razoumov 1998) gives an amplitude 4 or 5 times 
larger, which is still not enough to produce the warp, nor 
is the predicted direction of the warp in agreement with 
the observations. 

2.3. Example of application: magnetic torque 

Another example can be calculated: the case of magnetic 
forces (Battaner et al. 1990; Battaner & Jimenez- Vicente 
1998), in which the force per unit volume which produces 
the torque for a dp = 45° inclination field is: 



B 2 sin(2(ftp - a{R)) 
16ttZ 



)k. 



(17) 



where L — 1 kpc is the adopted value for the characteristic 
length in which galactic regions dominated by the galac- 
tic magnetic field become dominated by the extragalactic 




Fig. 3. Milky Way warp maximum height as a function 
of the radius due to gravitational tidal effects of the 
Magellanic Clouds (solid line). Dashed line stands for 
observational data of the northern warp (Burton 1988) 
rescaled to i?o = 7.9 kpc. 

magnetic field. Hence, the torque over a ring between R 
and R + dR is: 

rdR = 2h z R dR d<t>rxF (18) 
Jo 

B 2 R 2 h z sm(2(6 P -a(R)))dR / , . 
= (smcppi - cos (ppJli 

where h z = 0.1 kpc is the scale height of the disc. 

With the same Galactic model as the previous subsec- 
tion, we obtain that the amplitude of the magnetic field 
is 

B-1.4^G, (19) 

in order to obtain a coincidence with observational Burton 
(1988) data. The precession angular velocity is ui p = 5.1 x 
1CT 18 rad/s. 

There is a good agreement for the curves in Fig. |], 
which was noted by Battaner & Jimenez- Vicente (1998). 
The intensity of the field required to produce the warp is 
perhaps somewhat high, but we will not discuss in this pa- 
per the possibility of the existence of such a field. Kronberg 
(1994), for instance, argues that the value of of the inter- 
galactic magnetic field can be as high as this. We remark 
only that the possibility of warps generated by extragalac- 
tic magnetic fields should not be taken lightly, although 
the high value on the required field is perhaps question- 
able. This value is in agreement with that of Battaner & 
Jimenez- Vicente (1998) although their calculation method 
is much simpler. Binney (1991) obtained a higher required 
value for B, with a difference of an order of magnitude, in 
part due to his adoption of an expression for the magnetic 
force different from expression ( |l7| ) equivalent to a much 
lower value of L. He also adopted an insufficiently precise 



Fig. 4. Milky Way warp maximum height as a function 
of the radius due to an extragalactic magnetic field of 
intensity 1.4 /iG (solid line). Dashed line stands for obser- 
vational data of the northern warp (Burton 1988) rescaled 
to i?o = 7.9 kpc. 

approximation of an axisymmetric potential not distorted 
by the warp. 

3. Generation of warps by an intergalactic flow 

3.1. Accretion of matter onto the galactic disc and 
warps 

The idea we are suggesting here is not totally new. It is 
based on ideas about infall of intergalactic matter previ- 
ously explored by other authors (Kahn & Woltjer 1959; 
Binney & May 1986; Ostriker & Binney 1989; Jiang & 
Binney 1999). Kahn & Woltjer (1959) first suggested that 
an intergalactic matter flow could bring about the warp; 
however, their rough calculations in fact derive a pressure 
gradient transmitted to the disc by means of a hypothet- 
ical halo compressed by a subsonic massive wind and do 
not specify the mechanism of generation of the warp. The 
representation they use is very simple, and quite differ- 
ent from that presented here. The dominant response to 
the intergalactic medium ram pressure of a disc moving 
through it would be axisymmetric, taking the form of a 
rim rather than a warp (Binney 1992), although Kahn & 
Woltjer's estimates of the amplitude in fact agree with 
our values and their work was a first indication of a pos- 
sible mechanism to explain warps. Ostriker & Binney's 
(1989) suggestion of how cosmic infall can bring about 
warped galactic discs is rather qualitative and they did 
not proceed to a quantitative analysis. Subsequent papers 
(Binney et al. 1998; Jiang & Binney 1999) use a model 
based on the infall onto the Galactic halo and do not con- 
sider the direct infall onto the disc. The latter is the nov- 
elty here: we take the idea of cosmic infall, and explore 
the torque it creates when it collides transmitting its an- 
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gular momentum, but we consider directly the interaction 
with the disc. We do not invoke the idea of a massive halo 
modulating the dynamical effects on the disc. 

The idea of cosmic infall has been considered in the 
context of CDM theory. Ryden & Gunn (1987) and Ryden 
(1988) have shown that half of the total angular momen- 
tum of any galaxy was contributed by material that fell 
in over the last third of a Hubble time. Galaxy formation 
theories require this infall. Furthermore, many observa- 
tions imply that there must be an infall of material to 
galaxies (Binney 2000): Local Group members approach 
each other, the high velocity clouds around the Galaxy 
have on the average a net negative velocity (Blitz et al. 
1999; Braun & Burton 1999), and others. 

There are good reasons to believe that the infalling 
baryonic matter is accreted directly by the disc rather 
than the halo (non-baryonic matter can escape completely 
or be captured by the halo). The principal supporting ar- 
guments here are those based on the observed chemical 
evolution (Ostriker & Binney 1989; Lopez-Corredoira et 
al. 1999). Significant accretion of metal poor gas is nec- 
essary to justify the observations concerning star forma- 
tion and metallicity distribution in the galactic disc, of- 
ten termed the G-dwarf problem (Tinsley 1980) and the 
details of time-dependent evolution of individual metals 
(Casuso & Beckman 1997; 2000). Recent results implying 
that this accretion has been constant, or even increased, 
during the disc lifetime, are found in Rocha-Pinto et al. 
(2000). Moreover, it is clear that a halo should trap ac- 
creted matter with low efficiency, since its mean baryonic 
density is very low. 

In a more general context, it should not be thought 
that the general secular infall of matter is the only fac- 
tor. Any cloud in the intergalactic medium whose orbit 
intersects the galaxy and is accreted by the disc provides 
a torque due to the interchange of its angular momentum. 
For instance, the exchange of matter between two galaxies 
can supply accretable intergalactic matter: an intergalac- 
tic flow is produced between the two galaxies. The HVCs 
(High Velocity Clouds) have been suggested (Blitz et al. 
1999; Lopez-Corredoira et al. 1999; Wakker et al. 1999a; 
Binney 2000) to be observable evidence of the material 
which is continuously falling onto the Galactic disc. 

3.2. Description of the flow 

An intergalactic flow can be described as a beam of parti- 
cles which comes from infinite distance towards the galac- 
tic disc with velocity Vq. Each particle of the beam follows 
a trajectory which is not a straight line, due to the gravita- 
tional attraction of the galaxy, until it reaches the galactic 
plane (z — 0). As it intersects the plane, it collides with 
the gas of the galaxy and remains trapped in the disc. A 
torque results from the angular momentum contribution 
of this particle to the disc. The net torque over a ring of 
the disc of radius between R and R + dR will be the sum 
of the torques produced by all the particles of the beam 



which collide with the galactic disc at distance between R 
and R + dR from its centre. 

The total angular momentum with respect to the cen- 
tre of the galaxy which is transported by a cylindrical 
beam of these characteristics with axis crossing the cen- 
tre of the galaxy is zero. Is the total angular momentum 
deposited in each ring also zero? The answer to this ques- 
tion is "no" , and this is the key to warp generation by 
an intergalactic flow. The net momentum transferred to 
each ring is non-zero because, for a general case where the 
net flow is not perpendicular to the plane of the galaxy 
nor isotropic, the particles which fall onto a given ring do 
not come from a single cylindrical shell of the flow. The 
particles are redistributed, and so is their angular momen- 
tum, due to the gravitational interaction with the galaxy. 
Therefore, the impact parameter of each infalling particle 
is not the same for all the particles which cut a given ring 
(if it were the same the net contribution of the angular mo- 
mentum would indeed be zero); this will be analytically 
expressed in the equations given below in the present sec- 
tion. In particular, the variation of the impact parameter 
for a ring as a function of azimuthal angle is expressed in 
eq. ©. 

The calculations in subsections below are, perhaps, a 
little complicated to follow. However, the qualitative de- 
scription of the physical system whose variables we will 
calculate is not difficult to understand. Figure || gives a 
pictorial description which should be helpful. There is a 
net torque because the set of particles which fall into the 
ring between radii R and R+dR comes from a non-circular 
ring (which is not in fact elliptical). The transformation 
of this non-circular ring to the circular ring in the disc 
plane is effected by gravity. The velocity of the impacting 
particles varies with azimuth in the galactic plane ring. 
We wish to calculate the total angular momentum which 
is transported by matter in the non-circular ring which 
falls into the circular ring. Conservation of angular mo- 
mentum implies that the angular momentum transmitted 
to the circular ring must be the same as that in the non- 
circular ring. Nevertheless, calculation of this angular mo- 
mentum is not easy because the geometrical shape of the 
non-circular ring is not easy to describe analyticaly. This 
gives rise to the rather tedious calculations in subsections 
below. 

Once a torque is produced over each ring a warp is 
generated to compensate the differential precession of the 
rings, as we have seen in the previous section. We have ex- 
plained the mechanism of generation of warps when any 
kind of torque due to external forces is applied over the 
disc, in §|^. Now, in this section, we carry out the calcula- 
tion of the torque due to the intergalactic flow. 

3.3. Torque due to collision with a particle 

In Fig. ^, we give a graphical representation of the galactic 
disc centred at O, and the trajectory of a particle which 
intersects the disc. A particle comes from infinite distance 




vq = i>o(cos0o cos0o* + sin^ocosfloj + sinflofe)) (20) 



Fig. 5. Graphical representation of the infall of an inter- 
galactic flow (non-circular ring, above) which produces a 
net torque on a circular ring in the galactic disc (below). 



and it follows a trajectory which crosses the Galactic disc 
in some point Q whose distance from the centre is R 
and angle with respect to the x-axis (defined as the line 
"Galactic centre-Sun" with x negative towards the Sun) 
is d>, i.e. 



with velocity v a in the normalized direction e given by qq = r = r cos i + r q s i n( f) j 



(21) 




Fig. 6. Graphical representation of the Galactic disc in the 
coordinates system xyz, and the hyperbolical trajectory of 
an infall particle. 



A minor order correction for small warp amplitudes, 
that we also take into account, is the variation of angle of 
the flow in the warped rings. We must bear in mind that 



h(R)=9 (R = 0)-a(R). 



(22) 



The same trajectory is represented in Fig. [7] in the 
plane of the orbit x'y' . It is assumed that the trajectory 
is a hyperbola typical of a two-body gravitating system 
where the heavier body is the galaxy whose mass is M ga i 
concentrated at the point O. Some minor effects due to the 
dispersion of the mass throughout the disc are expected 
but they are negligible if R is larger than several disc scale 
lengths (i.e. greater than R « 10 kpc). The orbit is a 
hyperbola because the energy of the system is positive, 
since the velocity at infinite distance, |t>o|> is greater than 
zero. Therefore, the plane of the orbit will be determined 
by the independent vectors r and Vo and the equation 
of the orbit in the x'y' system is (r and (3 are the polar 
coordinates in this system; see Fig. ^): 



e 

— — = 1 + ecosp, 
A r 



(23) 



where e is the eccentricity of the orbit, 




G M, 



gal 



(24) 
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and 



A = 



G NL 



al 



1 



i 



(25) 



GM ga i e 

The impact parameter is b and the net asymptotic an- 
gular deviation 7 (see Fig. |^) is given by 

GM gal 



7 

tan- 



H 2 



(26) 



The determination of the point of intersection of the 
orbit with the disc of the galaxy is a simple trigonometric 
problem. From a triangle shown Fig. |?], we can derive: 



7T 7 

PQ = 2 + 2 



e Q = cos(v ,tq) 



l (eoQ), 



cos(6* ) cos((/>o — <f>). 



(27) 
(28) 



From these expressions, together with fl23|), (|24j) , (|25j) 
and ( p6| ) , the radial galactocentric distance R of the point 
of orbit intersection with the Galactic plane is: 



R = r Q 



b 2 vl 



bvl 



1 e 0Q 



GM gal (l-e 0Q ) 



The angular momentum of the particle with mass dm 
is constant along its trajectory, 



dm vob- 



Tq X Vq 



(30) 



r Q v \sm(v o ,r Q 

and the torque produced by a particle which transmits its 
angular momentum to the disc is: 

Vq x v dm 
~ ~d~t 



dJ dm 

-dT = 



rQ V 



1 



= ^- Vo b(l-el Q r^ 



-0Q 



x [sin (f> sin 9qx — cos <f> sin 6*0.7 + cos $0 sin(</>o — <f>)k] . (31) 

The new material is stopped by the friction with the 
disc and its angular momentum added to the ring. Note 
that although the ring is not really a rigid body it behaves 
like one. A single particle in orbit around the Galactic 
centre can resemble the dynamics of the rigid body; it 
carries the increase of angular momentum and the orbit is 
distorted according to the applied torque. 

The flow is stopped by the friction with the gas, so 
the angular momentum is, at first, transmitted to the 
gas. This does not mean that gas disc warps while stellar 
disc does not. Indeed, it is expected that stars in the ring 
feel the gravitational torque due to the gas rings and are 
dragged by them. There may be some difference between 
the stellar warp and the gas warp due to this lag in the 
dynamics, but the difference is likely to be small. If the 
stellar disc were demonstrated to be less warped than the 
gas warp, it would be evidence in favour of either this 
theory or the theory of the intergalactic magnetic field as 
the generator of the warp (Porcel et al. 1997), or indeed 
of any theory in which the external torque affects directly 
the gas disc rather than the stellar disc. 




Fig. 8. How the polar angles (j>b in the original flow and 
(29) </> in the plane of the galaxy are related in terms of a 
spherical triangle. 



3.4. Total torque due to collision with a particle beam 

Expression ( |3l| ) gives the torque produced by a particle 
which falls to the disc with an impact parameter b and in- 
tersects the disc at tq . If we want to know the total torque 
produced by all the particles which come with any b and 
intersect the disc at a distance between R and R+dR with 
any azimuth <p, we have to integrate over all the particles 
of the beam which fall in this ring. The whole beam is 
then represented by the varying in the plane perpendicu- 
lar to Vq the initial (at infinite distance) position of the 
falling particle, whose polar coordinates are b and 4>b (see 
Fig. |8j). Thus, the total torque exerted over the ring with 
radius between R and R + dR is 



T(R)dR = [ d<t> b [ 
Jo Jo 



dbb 



;R<Rq <R+dR 



,dm 



x fjjT v o b (! ~ eog(0) 2 ) 1/2 [sin0sin6) o i 



cos sin #0.7 -I- cos 0o sin(0o — </>)&]], 



dm = pb vodt, 



(32) 
(33) 



where ph is the density of baryonic matter in the particle 
beam, assumed to be independent of b and fa. Any non- 
baryonic matter in the inflow would not be trapped in 
the disc so it should not be taken into account within 
the total mass of the flow for the purpose of computing 
the torque. In the notation, r stands for the torque per 
unit galactocentric radial length. Note that fit, is the polar 
angle in the plane perpendicular to v and is different from 
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the polar angle 4> m the Galactic disc. The relationship 
between the two is: 



cot( 



>b - Pbn 



COt(</> - 4> n ) 



I sin O 



(34) 



according to the general formula relating the angles in 
the spherical triangle of Fig. ||. <p n and 0b„ are the polar 
angles of the galactic disc and the plane perpendicular to 
Vq respectively corresponding to the node where the two 
planes intersect. We can choose for convenience the origin 
of the angles (f>b such that: 

4>bn = 0, (35) 

and the angle of the line of nodes in the galactic disc is: 

cf> n = 0o ± tt/2. (36) 

We change the variables of integration in the expres- 
sion (52]) to R and <f> (the Jacobian of the transformation is 
l^^r 1 ! we neglect f| due to the dependence 6 Q (R)) 
and obtain 



t(R) 



^(l-epQ^) 2 )- 1 / 2 
I sin O I Jo 1 + sin 2 (0 O - 4>) (sin" 2 9 - 1) 

x [sin sin Ooi — cos sin 80 j + cos 60 sin(0o — </>)&] 

<%(#, 



Si? 



-b(R 



2. 

5 V/ 5 



where b(R, </>) is derived from ( ^9|) 
6(i?,0)-ii?^l- e 2 Q (0) 



(37) 



(38) 



+ ^ i i? 2 (l - e 2 Q (0)) + R G M gal v 2 (l - e QQ m. 

If we make a new change of variables in the inte- 
gral, substituting x = eog(</>) for <j) according to (f28|), 
and simplify the expression (some terms cancel because 
of the antisymmetry in the interval (0 — (f> ) £ (0, it) and 
{(f) — 4>o) e (7r, 2tt)), we obtain 



2p b ?;o sinflo , . 



sm Co cos^ 6*0 



( cos9 ° db(R,x) Lfn a 

x / a; — b(R^x) 

a OR 

- COS fg 



(39) 



VI - s 2 v/l - a; 2 / cos 2 O [1 + (1 - a; 2 / cos 2 6) )(sin" 2 () - 1)] 

The integral is positive. The direction of the torque 
is ±(sin0o* — cos 00.7'); the sign is '+' when 80 < and 
'-' when 80 > 0. That is, the torque is in the disc and 
perpendicular to Vq. 

The direction of the infall of a particle flow is the 
same as that produced by the gravitational torque when 
e = ep. although the amplitude is rather different. Note, 
for instance, that the amplitude in (p9|) does not depend 



on the disc density. Therefore, the effect produced by the 
infall of the particle beam is similar to that due to gravi- 
tational effects. 

The torques between rings do depend on disc densi- 
ties (expression (H) is proportional to a(R)a(S)). This ex- 
plains why the disc will warp significantly only where its 
surface density is low (Ostriker & Binney 1989) which, in 
practice, means at its outer edge. For the inner disc, a is 
high enough to provide strong torques between two rings 
with a small angle a^s of separation which can compen- 
sate the difference of precession with respect to the av- 
erage. However, external rings must separate further to 
compensate these differences. 

Note that there is no z-component of the torque. 
Taking into account that the total mass of the ring is 
increased, this means that the angular velocity of rotation 
should decrease and, therefore the radius of the orbit will 
be reduced (the increasing mass of the inner galaxy will 
also tend to reduce the radius) and the galaxy will con- 
centrate further material in the inner regions. We will not 
study these aspects further in the present paper. 

The limit of low initial velocity (^r^f 1 « 1) 
gives a proportionality 



I 4G JV/ gal j 



lim T {R) Ph 



(GM gal )3/2 



(40) 



x (sin 0oi — cos (f>oj). 

This means that a low velocity beam provide a stronger 
torque. This effect can be seen due to the increased cur- 
vature of the hyperbolic trajectories for low velocities. A 
ring will then accrete flow particles over a wider range of 
impact parameters, if the velocity of the beam is lower. In 
the extreme case of vo = all the particles fall to the cen- 
tre, so R — for all cases with a finite impact parameter 
and there is no divergence (i? cx Vq from (f29|)). Although 
we have integrated b, the impact parameter, from zero 
to infinity, in a real case b is limited to a finite value. 
Expression (40) means that at lower flow velocities the 
disc accretes particles from a greater fraction of the total 
stream of particles, and this is the reason why extra angu- 
lar momentum is deposited. Again, as emphasized above, 
we must note the angular momentum deposited within a 
given annulus is due to accretion of particles from a non- 
axisymmetric volume of space whose integrated angular 
momentum is thus non-zero. 



3.5. Mass accreted by the galaxy 

The total accretion rate due to this infall is 
dM r27r 



dt 



A h ,; h " 

db b— = p h v 

JQ; i?<fl max at J 



R ^dR b{RA" dh[R 



dR 
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-PbV 



sin O 1 cos 0, 



o J- 



dxb{R n 



^Jl-x 2 /cos 2 O [1 + (1 - x 2 / cos 2 O ) (sin" 2 6> — 1)] 



.(41) 



^1 -x 2 /cos 2 9 Q [1 + (1 - x 2 / cos 2 (9 )(sin- 2 9 - 1)] ' 

This acceleration is compensated by the vertical grav- 
itational acceleration, so: 



This rate is positive since the integral is negative 
(b{R,-x) > b(R,x)). 

3.6. Transfer of linear momentum 

It is important at this stage to consider the transfer of lin- 
ear momentum. While the transfer of angular momentum 
produces a integral sign shape (component m = 1 of the 
galactic warp), as shown above, the linear momentum pro- 
duces a cup-shaped deformation of the disk (component 
m = of the galactic warp). The numbers in the section 
§[| will show that this effect is indeed quite small, and the 
m = 1 component is predominant. 

The points of equilibrium between the vertical gravita- 
tional force and the vertical force due to transfer of linear 
momentum give the deformation of the disk. The verti- 
cal linear momentum due to the accretion of a particle of 
mass dm is: 



Pz 



dm v sin(i?Q ,v)yl — cos 2 9q sin 



(42) 



The vector v is the velocity at the impact point of the 
disc (Rq). The last two factors account for the projec- 
tion of the velocity onto the vertical axis. It is not easy to 
visualize the origin of these factors; but they can be un- 
derstood by reference to Figures || and 0. The radial and 
azimuthal linear momentum would produce some distor- 
tion of the orbits within the disc, which is not the subject 
of the present paper. 

From the conservation of the angular momentum, we 
have: 



\J\ = dm vob = dm v Rq sin(i?Q, v) 



(43) 



From expresions ([42]) and (|43|), the vertical force of a 
particle of mass dm is: 

"-WW.-*. («) 

The total vertical force on all the particles of the beam 
which fall into the ring of radii R and R + dR is (from now 
on, F z stands for the vertical force per unit galactocentric 
radial length): 



F z {R)dR= I acp b 

10 JO; R<Rq <R+dR 

— cos 2 9q sin 2 (0o — 4>) 

/•cos 9q 



M h ^R^ (45) 



= 2pw4dR r° se ° c , db(R,x)b(R,xf 
| sin O I cos O 7- cos e dR R 



F z (R)dR GM gal (i?) z 



2ira(R)RdR R 2 R' 



(46) 



which implies, for a mass derived from the rotation curve, 
that 



F Z (R)R 



2Trv mt (R) 2 a(R) 



(47) 



This monopolar approximation would be exact if the 
distribution of mass were spheroidal or elliptical. For the 
disc, whose contribution is dominant in the torque, the 
contribution of the gravitational force is different: 

F (R) /"flmax 

- — 2 -. —. — — Gz I ds s cr(s) 

2tt jj, 



o [s 2 + i? 2 -2si?cos(/> + z 2 ] 3 / 2 



Gz ds s ens) / -f-= — 5 — — rm^- 

Jo J Q [s 2 + R 2 - 2sR cos cj)] 3 / 2 



(48) 



The last approximation is for z small, compared to (R—s); 
this is a good approach for large R. This is also an approx- 
imation in another sense: it does not take into account the 
distortion of the disc produced by both m = and m = 1 
components of the warp. In any case, most of the mass is 
in the inner rings which are not distorted and this effect is 
small. Hence, z is inversely proportional to the attraction 
of the disc even for a non-monopolar approximation. 

The disc contribution is something larger than its 
monopolar contribution, so z is lower than (f47|). This 
means that the expression (^) gives, approximately, a 
maximum limit for the distortion of the disc, i.e. 



z < 



F Z (R)R 



2TTv mt (R) 2 a(R)' 



(49) 



This approximation is good enough when the effects of 
the transfer of linear momentum are nearly independent 
of those due to the transmission of angular momentum. A 
more accurate model would refine these approximations 
and calculate both the m — and m = 1 distortions 
simultaneously. We will see later that the effect of the 
transfer of linear momentum is small compared with the 
transfer of angular momentum and, therefore, we do not 
need to make these more refined calculations for this effect, 
but only for the integral sign warp, which is dominant. 

Is it counterintuitive to have m = 1 deformation larger 
than m = deformation? We can argue that it is not. 
Imagine a galaxy in which practically all the mass is in the 
centre, so that the disc has a negligible mass. How big is 
the countertorque of the galaxy to compensate the exter- 
nal torque? The answer is: zero (or nearly zero), because 
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only the dipole and higher order components of the gravi- 
tational attraction produce a countertorque. The monopo- 
lar component does not create countertorque. In this case, 
what will be the equilibrium condition in which the torque 
is compensated by the countertorques of the disc? It is an 
"m = 1 deformation of the disc" (i.e. a warp) whose am- 
plitude tends to infinity (as the mass of the disc tends to 
zero). In fact the limiting case will not be at infinity, since 
we are talking about an angular amplitude. The ring will 
take the angle of incidence of the wind as its maximum 
limit. So what happens to the "m = deformation"? This 
deformation is nothing like as strongly dependent on the 
distribution of the mass and there is a net finite counter- 
force of the galaxy to compensate the external force; the 
amplitude of the m — cup-shaped distortion is thus fi- 
nite. Therefore, we have a very high amplitude for m = 1 
and the amplitude of m = is much smaller. If the density 
of the wind tends to zero, the amplitude of m — distor- 
tion tends to zero, but the amplitude of the m = 1 dis- 
torsion tends to take up the incidence angle of the wind. 
We conclude that this is a very clear case in which the 
m = 1 distortion amplitude is much bigger than the m = 
amplitude and, therefore, it is possible to have a m = 1 
deformation larger than the m = deformation. 

3. 7. Direction of the warp 

One relevant aspect of the present hypothesis is that it al- 
lows us to relate the direction of the warp to the direction 
of the inflow, Vo, with respect to the centre of the galaxy. 

The torque between two rings, dTs(R)dR dS, is pro- 
portional to (sin 4>oi — cos <froj). This can be demonstrated 
by use of the equation (|^), or by a simple analogy with 
the gravitational torque of one particle (§2.2): the second 
ring is the equivalent to an integration of the particle po- 
sition over the second ring azimuthal angle. Likewise, the 
external torque is proportional to (sin^o* — cos0oj')- This 
means that the warp direction is the same of the projec- 
tion of Vo over the disc, and the precession is around the 
axis parallel to Vo. 

The amplitude of the collisional torque, (p9|), is larger 
for larger R, since b 2 increases with R, so the outer rings 
have an excess angular velocity of precession compared to 
the average. However, the outer rings receive a torque from 
the inner rings, which is proportional to the collisional 
torque when (f>p — cf>o, where ep is the position of the 
maximum height above the plane of both inner and outer 
rings (i.e. the direction of the warp), but must act in the 
opposite direction to counteract the above excess. Hence, 
the outer rings should be oriented towards positive z if 9q 
is positive or towards negative z if 0$ is negative. 

From the above considerations, we deduce that the 
warp is oriented in the direction parallel to the projection 
of erj in the disc. The azimuthal angles of the maximum 
and minimum heights (z) of the warp are </>o and fa + tt 
respectively if 8o > 0. If 9 < 0, the maximum height is 
at azimuth (f>o + n while the minimum is at 4>q . Figure ^| 




O 



Fig. 9. Graphical representation of the orientation of the 
warp with respect to the velocity of incoming flow. The 
azimuthal angle of the warp, which is not plotted, is the 
same as the azimuthal angle of Do, i.e. (f>o. 

shows a graphical representation of these orientations in 
the plane which contains the vector Vq perpendicular to 
the disc. Note that the orientation of the warp does not 
change if the velocity is — Vq instead of Vq. 

4. Warp in a typical spiral galaxy due to 
intergalactic accretion flows 

The proposed hypothesis of formation of warps can be 
compared with our observational knowledge about galax- 
ies and the intergalactic medium. We must bear in mind 
that the calculations developed in this paper are based on 
fairly crude approximations. We have considered particle 
trajectories as hyperbolic, neglecting the gravitational ef- 
fects of the extended disc, we have considered a collimated 
infinite beam, constant density of the intergalactic flow, 
etc. Therefore, we should not expect the theory to repro- 
duce exactly all the fine details of a warp. However, we will 
show that our model reproduces quite well the observed 
warps using a standard model of the disc of a typical spi- 
ral galaxy (for instance, the Milky Way) and with entirely 
reasonable parameters for the intergalactic flow. Our re- 
stricted goal here is to offer a general model for warps (a 
warp amplitude a(R) which is close to zero out to some 
specific R, but rises rapidly at higher increasing i?), and 
to test it using the parameters of a typical spiral galaxy 
such as the Milky Way. 

4.1. Amplitude of the warp 



The solution of ( |12| ) by means of the numerical method 
explained in appendix |A| gives a(R) for the galactic warp. 
It will depend on the adopted galactic model. Our pur- 
pose is to demonstrate that the order of magnitude of 
the amplitude of the warp is approximately the observed 
amplitude, so we use only one model. Varying the param- 
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eters of the disc will, of course, lead to variable results 
but these would not be quantitatively very different. The 
surface density of the disc we adopt is given by @ and 
the rotation velocity is given by ([l6|), both corresponding 
to the Milky Way. 

In the linear regime, for a low amplitude of the warp 
and low velocity vq, in the limits (EG) and (o), the pro- 

M 3/2 U 

portionality followed is a(R) oc lo p oc — gal , in which the 
dependence of a{R) on R is not specified. Specifically, for 
the adopted disc model and \6q\ = 45°, we obtain in this 
linear regime that the angle of the warp is (0& = 0.02 h~ 2 ; 
Schramm & Turner 1998) 



a(2i? ) ~ 5.2 x 10~ 34 



p b (kg/m 3 )M gal (kg) 3 / 2 



Data from Burton (1988) 

Prediction from the infall model (m=0) 

Prediction from the infall model (m=1) 
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Burton (1988) gives an observational value of 
a(2Ro) w 0.14 rad in our Galaxy which, if we take the 
value at 2R$ to normalize the amplitude of the warp, leads 
to 



Fig. 10. Milky Way warp maximum height as a func- 
tion of the radius in the hypothesis of continous accre- 
tion of intergalactic matter with pb = 1-1 X 10 -25 kg/m 3 , 
|f?o| = 45° (m = 1, solid line). For comparison, a maxi- 
mum limit of the effects of the linear momentum incre- 
rad/s.(51)nent over the rings are shown (m = 0, dot-dashed line). 

Dashed line stands for observational data of the northern 
warp (Burton 1988) rescaled to Rq = 7.9 kpc. 
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6.2 x 10~ 18 rad/s = 1 cycle in 32 Gyr. 



(52) 



(53) 



The precession is far too slow to be significant, and much 
slower than the rotational velocity of the Galaxy. 

We adopt a total mass of the Galaxy of M ga i = 2 x 10 11 
M Q (Honma & Sofue 1996) to calculate the curvature of 
the hyperbolic trajectories. This estimate includes bulge, 
disc and spheroidal component masses. Hence, 



Pb ~ 290 fibPcrit 



100 Km/s 



(54) 



The only free parameters with observational uncertain- 
ties are the mean density of the intergalactic inflow and its 
velocity which are related by the last equation. Assuming a 
typical velocity, in galactocentric coordinates, of vq ~ 100 
km/s (the limit of low velocity adopted in ( |40| ) is valid if 
R is much less than 300 or 400 kpc, which is certainly the 
case since we take R < 16 kpc), we find: 



p h ~ 290 ftbPcrit = 1-1 x 10~ 25 kg/m 3 , 
which is equivalent to 
n m ~6x 10~ 5 / cm" 3 , 



(55) 



(56) 



where / is the fraction of HI in the total baryonic mass 
of the wind. This value coincides with the most probable 



value estimated by Kahn & Woltjer (1959) for the mean 
density of intergalactic matter yielding dynamical stability 
for the Local Group of galaxies. Lopez-Corredoira et al. 
(1999) also calculated a total intergalactic mass in the 
Local Group around 2 x 10 12 M Q , which in a volume of 
~ 1 Mpc 3 gives a mean density around 10 -25 kg/m 3 . 

The shape of a(R) w \z\/R for these values would be 
that given in Fig. |l(j. The prediction of this model agrees 
with many of the features of the observed warp (Burton 
1988; 1992): it predicts a flat disc which is not significantly 
warped for values of R less than ~ 1.3 Ro — 10 kpc, and 
a warp whose amplitude increases rapidly at larger radii. 
The shape of the warp depends on the detailed radial vari- 
ation of a(R) and we have used a rough exponential es- 
timate of this dependence, which is really more complex 
than a simple exponential function with a constant scale 
length. In any case, it is not at this stage worth using a 
more realistic model since our approach is only an approx- 
imation, given the simplifications we have employed in the 
analytical calculation of the torque. The functional shape 
of a(R) resembles that of the observational data, within 
the values of R less than ~ 16 kpc from the centre. Non- 
linear effects have a significant influence at larger radii as 
well as departures from the simple ring model of the disc. 
Nevertheless, we should note that the dependence of x cxt 
on R does not affect too much the form of a(R) and other 
external torques may also produce similar shapes, as we 
have seen in Figures || and ^. 

The maximum dist ortio n of the disc due to the transfer 
of linear momentum (§ |3.6| ) is calculated using the expres- 
sion (fl9j). The results are also plotted in Fig. [To|. In this 
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Fig. 11. Milky Way warp maximum height as a function 
of the inclination of the incoming flow with respect to the 
galactic plane (9q) in the hypothesis of continous accretion 
of intergalactic matter: m = 1 (integral-sign warp due to 
a torque), solid line; m = (cup-shaped distortion due to 
a force), dashed line. 



case, z is always positive or always negative, a cup shape 
rather than an integral sign shape. From the numbers plot- 
ted in this figure, it can be concluded that this distortion is 
small compared to that produced by the torque, i.e. the in- 
tegral sign shape will be predominant in the galactic disc. 
We think the predominance of torque effects is due to the 
smaller mass concentrated in the disc in comparison to 
the total mass of the Galaxy, which is mostly responsible 
of the counter-torques which reduce the heightness of the 
warp m = 1, while the whole mass of the Galaxy produces 
the counter- forces which reduce the distortion m = 0. In 
some galaxies, a mixture of cup-shape and integral sign 
shape could be produced by this mechanism. This would 
be a possible explanation for the asymmetries between the 
two sides of the warp, even for our own Galaxy. It will de- 
pend, among other factors, on the angle 9q. For 6*o = 90°, 
the transfer of vertical linear momentum would be maxi- 
mum while the torque would be zero. Therefore, there is 
a probability, although small, of finding cup-shaped dis- 
tortions of galactic discs rather than integral-sign warps. 
In Fig. [ll], we show that the cup-shaped distortion is pre- 
dominant for \8q I >^ 85°, which translates to a proba- 
bility that a cup-shape is dominant of <^ 4 x 1CP 3 , for 
spiral galaxies like the Milky Way. It would be of inter- 
est to make a statistical search over a major sample of 
warped galaxies to see how many galaxies present appre- 
ciable cup-shaped distortion. The study of the fraction 
of galaxies which present irregulatities within the integral 
sign warp would also give some clue about the validity of 
this hypothesis. 

Near the limit 6$ = the amplitude of the warp goes to 
zero quickly. Although, the figure Fig. |ll| shows apparently 



a vertical tangent in the limit 9q = 0, that is not the case: 
there is a horizontal tangent, as the limit of (|3^) when 
9o is very small is proportional to 9q. However, this limit 
manifests for very small angles and this is the reason for 
the apparent vertical tangent in [H]. Indeed, the values 
are: \z\(9 = 10°) = 2.45 kpc, |z|(0 o = 5°) = 2.24 kpc, 
\z\{9 = 2°) = 1.79 kpc, |z|(0 o = 1°) = 1.24 kpc, \z\(9 = 
0.5°) = 0.59 kpc, |z|(0o = 0.25°) = 0.19 kpc, \z\(9 = 
0.1°) = 0.034 kpc,..., z{9 = 0) = kpc. 

We can add some further comments about the stel- 
lar warp. In the Milky Way, the OB stars follow the gas 
(Porcel & Battaner 1995). But it is not just the young 
population, recently formed from the warped gas, which 
shows a stellar warp, but the whole population of stars 
(Porcel et al. 1997; Dehnen 1998). The whole population 
of stars of the Milky Way projected onto the sky appears 
less deviated from the plane than the projected gas , but 
this does not necessarily represent a difference between 
amplitudes of the gaseous and stellar warps. Rather, it 
may well be due to the cut-off of the stellar population at 
a smaller Galactocentric radius than the gas, around 15 
kpc (Porcel et al. 1997). The stellar disc is clearly warped, 
perhaps somewhat less than the gaseous disc although not 
obviously so. As pointed out in §3.3, if the stellar disc 



were demonstrated to be less warped than the gas warp, 
it would be evidence supporting either this theory or the 
theory of the intergalactic magnetic field as the generator 
of the warp (Porcel et al. 1997). In fact, this would be pre- 
dicted by a theory in which the external torque directly 
affects the gas disc and only indirectly the stellar disc. 

4.2. Possible scenarios 

Possible scenarios for the intergalactic flow postulated 
here are: 

1. The galaxy is passing through a continuous intergalac- 
tic medium, i.e. the velocity of the flow is due to the 
relative motion of the galaxy with respect to the rest 
frame of the intergalactic medium. In this case, the 
flow may be well approximated by a beam of infinite 
extent. This may be a representative scenario for most 
galaxies, and could explain why most spiral galaxies 
are warped. Asymmetries in the warp could also be ex- 
plained in terms of a transfer of both linear and angular 
momentum, although the cases when cup-shaped dis- 
tortions dominate has low probability. We think that 
this first item is the most plausible explanation. 

2. The galaxy has an interacting companion and some 
exchange of material is produced due to tidal effects. 
The stream of material goes from the companion to 
the main galaxy, and this flow could also generate a 
warp. The direction of the wind would not be constant, 
because the companion is orbiting around the main 
galaxy so the warp would not be steady. This scenario 
would give asymmetries in the warp, apart from that 
coming from the m = component, which would be 
due to the non-symmetric form of the infalling mate- 



16 



Lopez-Corredoira et al.: Warps & intergalactic flows 



rial, with finite impact parameter, or clouds which in- 
tersect the centre of the galaxy in a non-axisymmetric 
way. Reshetnikov & Combes (1998) find some corre- 
lation between the frequency of warps and the inter- 
action with other galaxies. An accretion inflow due to 
the exchange of material with a companion might ap- 
pear to be undermined by the observations of warped 
galaxies which are apparently isolated. However, this 
phenomenon is itself being challenged by the discovery 
of some companions to galaxies which had been con- 
sidered as classical examples of isolated galaxies (for 
instance, NGC 5907; Shang et al. 1998). 

While the present paper was in the refereeing process, 
a work by Garci'a-Ruiz (2001) has been published. This 
is an interesting work which analyzes in detail 26 edge- 
on galaxies in radio and optical. Garcia-Ruiz (2001) finds 
that 20 galaxies are warped, two of them present U-warp, 
and 7 with a warp of only one side. This asymmetric cases 
can be explained by the present theory as a combination 
of m = and m = 1 distortion and it is interesting to 
note that up to now there is not any alternative explana- 
tion for the U-warps as well as the asymmetric warps. He 
also find that the frequency of warps and its amplitude is 
dependent on environment. It is even more interesting to 
note that the most isolated galaxies are more frequently 
warped (although with less amplitude). "While this sug- 
gests that tidal interaction plays a role in warping, it seems 
likely that there are other effects at work that cause even 
quite isolated galaxies to warp." (Garci'a-Ruiz 2001). It 
seems clear that warping is due to something related with 
the environment rather that the intrinsic propierties of 
the galaxies, and something which is not related with the 
proximity of other galaxies. The accretion of intergalactic 
matter onto the disc seems a good candidate to explain 
these observational facts. 

The density of the intergalactic medium is very low on 
average, making it difficult to detect. However, it could 
well be that the High Velocity Clouds (HVCs; see the 
reviews in Wakker & van Woerden 1997; Wakker et al. 
1999b) are part or all of this material falling towards the 
Galactic disc from all directions. If these do produce a net 
torque it must be because there is a net galactocentric av- 
erage velocity of the complete set of HVCs, including the 
Magellanic Stream, the complexes and the HVCs associ- 
ated with more distant clouds infalling towards the Local 
Group barycenter (Blitz et al. 1999; Braun & Burton 1999; 
Lopez-Corredoira et al. 1999). The tidal Stream from the 
Saggitarius dwarf galaxy (Ibata et al. 2001) might also 
take part. An average hydrogen density for the flux of 
Uhi ~6x 10 -5 / cm~ 3 is in good agreement with the av- 
erage density of an HVC, (uhi) ~ 10 _4 -10 _1 cm -3 (Blitz 
et al. 1999; Wakker et al. 1999a), including the value of / 
which may be as low as 0.02 (Braun & Burton 2000) if the 
matter of the flow is baryonic. The mean density in the 
intergalactic clouds may be enough to create the Galactic 
warp. If the density of infalling matter were equal to that 
of an individual HVC, we would have a higher density than 



that required, by several orders of magnitude. However, we 
know that the intracluster medium is not filled by HVCs; 
these represent a low volume high density fraction; averag- 
ing over the complete medium with plausible filling factor 
can yield a net hydrogen density around rim ~6x 10~ 5 / 
cm~ 3 . 

The degree of dumpiness of the intergalactic medium 
is not well known. However, it is not of importance 
whether the flow is continuous or discretized in clouds. 
The warp will be produced by the average infall. Short- 
term fluctuations of the infalling density do not apprecia- 
bly affect the warp since the forces responsible for it have 
a very low amplitude and require a long time to produce 
or distort the warp. 

Using expression (|4l|), we can derive a total accretion 
rate of the Galactic disc out to i? max = 15 kpc of ~ 1 
M Q /yr for this density. This turns out to be of the order 
of the accretion rate required to resolve the G-dwarf prob- 
lem in our Galaxy as well as explaining a number of phe- 
nomena of chemical evolution which require the long-term 
infall of low metallicity gas (Lopez-Corredoira et al. 1999; 
Wakker et al. 1999a). The infall of 1 solar mass per year is 
enough to produce the warp because: 1) the external disc 
has a very low density and small forces produce consider- 
able accelerations; 2) the acceleration may in fact be very 
small in amplitude, but the period of time to produce the 
warp is large enough (order of Gyr), so this gives time to 
distort the galactic disc (in 1 Gyr, 10 9 solar masses are 
accreted which is a significant quantity of accreted mass) . 

This hypothesis could, in a general way, explain the 
possible alignments of the different warps of neighbour- 
ing galaxies (Battaner et al. 1991) if the flow velocity is 
similar around such galaxies within the same zone of the 
intergalactic medium. 

Whatever the structure and composition of the inter- 
galactic medium, it is clear that intergalactic space is by 
no means empty, and the accretion of this material by 
galaxies is likely to have been continuous during their life- 
times. The effects of this accretion can be detected in their 
chemical evolution as well as in their structure, as pointed 
out above. 

5. The effects of a very massive halo 

Although we already discussed the effects of the halo, we 
are going to clarify the use of the halo here as well as the 
effects which a very massive halo could produce. 

We may infer from the present calculations that we 
have found a mechanism which could explain both quali- 
tatively and quantitatively the generation of warps in nor- 
mal spiral galaxies, and this mechanism is the interaction 
between the disc and the infalling matter as well as the 
interaction of the disc with itself (rings interacting with 
other rings). The dark halo is included although it is not 
explicit in some calculations of this paper, but its effect is 
not very important for a rough calculation in which we are 
interested in the order of magnitude. In fact, our calcula- 
tions are not rough but almost exact, but the roughness 
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of the numbers obtained is due the uncertainties in the 
different parameters. 

Briefly, these are the effects which are treated here, and 
others which have not been treated but require treatment 
in future papers: 

— The dark halo is present in the total mass of the 
Galaxy: M ga i. The assumed value of the Milky Way 
mass of 2 x 10 11 M within a radius R ~ 20 kpc from 
the centre includes the halo, some of whose matter is 
not visible. The effects of the halo are also implicitly 
included in the calculation of M ga \(R) as a function of 
Vmt(R) from the rotation curves. 

— The interaction of the halo with the warp is not consid- 
ered, but we have made estimates to show that for this 
case the effect is small compared to the countertorquc 
of the disc. Therefore, we feel that the exclusion of halo 
effects in the countertorques is a fair approximation if 
we are interested only to find the order of magnitude 
of the warp amplitude. The countertorques of the halo 
are less important than those of the disc because the 
halo is much closer to sphericity. We showed in §|| that 
the contribution of the halo countertorque for R < 16 
kpc is less than ~ 40% of the disc contribution. The 
uncertainties in the halo mass distribution lead to an 
uncertainty similar to this value and we did not feel 
that much would be gained by adding the extra com- 
plexity. We know that the warp would be reduced by 
this effect, but this reduction is not the dominant term. 

— A very massive and very extended halo would intro- 
duce some variations in the numbers we derived al- 
though the mechanism would work qualitatively in a 
similar way. First, if it is very massive M ga \ would be 
larger and the amplitude of the warp would be corre- 
spondingly larger. We have shown that the amplitude 

3/2 

of the warp is proportional to M . This, however, 
was on the assumption that most of the mass is concen- 
trated within a radius less than ~ 20 kpc. If the halo is 
very extended, this approximation is not appropriate. 
We should then need to consider the mass distribu- 
tion of the halo instead of assuming a point mass to 
calculate the infall velocity of the clouds. The prob- 
lem would be much more complex because we would 
need a gravitational potential different from <f> oc 1/r 
and the trajectories would not be hyperbolae. Indeed, 
the mechanism we propose to form warps works with 
any potential, but the calculations are of course much 
easier for a oc 1/r law. Since we are interested in 
proposing a new mechanism and showing how it works, 
we think that these complexities should be left for a 
future paper. At present, we can say that a very mas- 
sive halo would increase the amplitude of the warp, not 

3/2 

proportionally to M ^ but some other power with ex- 
ponent less than 3/2. This is because the infall velocity 
of the clouds is increased as well as the curvature of 
the trajectories. It is important to note that a very ex- 
tended halo would not differ from a halo constrained 
within R < 20 kpc in the countertorques produced on 



the warp because only the mass in the ellipsoids in- 
ternal to the radius of the warp produce gravitational 
torque, assuming there is a constant ellipticity halo. To 
summarize, with a very massive halo our mechanism 
does operate, even in fact more effectively because it 
would require a lower intergalactic density to produce 
the same warp amplitude. 
— If the halo axis were displaced with respect the disc 
axis the disc would be pinched by the halo. This was 
the case studied by Ostriker & Binney (1989), Binney 
et al. (1998) or Jiang & Binney (1999) and is indeed 
a mechanism which will give rise to warps too. The 
scenario presented by these authors is different from 
the scenario proposed here. They assume the infall of 
the intergalactic matter onto the halo rather than the 
disc, and this produces motion of the halo with respect 
to the disc. We do not seek to challenge this in the 
present paper, rather to present a possible alternative, 
which might be complementary. The compatibility of 
the two mechanisms certainly merits effort to study 
further. There is no doubt that accretion into a halo 
which gives an offset in the rotation axis from the disc 
axis will give rise to a warp; if one accepted the possi- 
bility of a low density halo accreting intergalactic mat- 
ter, one should recognize the possibility of giving rise 
to warps by this mechanism. We have shown here how 
accretion directly onto the disc can also yield a warp, 
with parameters in the observed range. At the present 
stage of understanding the problem, we think that ei- 
ther or both mechanisms can act to produce warps if 
one accepted that the accretion of matter by the disc 
or the halo have the same plausibility. Our opinion is 
that our mechanism is preferable since the accretion of 
matter by the disc is more plausible than the accretion 
by the halo but this point is open for discusion. 



6. Conclusions 

We propose that galactic warps are produced by the reori- 
entation of the galactic disc structure in order to compen- 
sate the differential precession due to a torque generated 
by an external force. The external force might be the grav- 
itational interaction with a satellite, but in the case of the 
Milky Way, with the Magellanic Clouds as the satellites, 
there is not enough mass close enough to provide the ob- 
served amplitude in the warp. Magnetic forces could also 
produce the warp, but the intergalactic field would then 
be of the order of /iG. 

A simple model of an intergalactic accretion flow which 
intersects a galactic disc (or, equivalently, considering the 
galaxy as moving through the intergalactic medium) can 
explain the existence of warps in the galaxy if the mean 
density of baryonic matter in the medium is around 10 -25 
kg/m 3 and the infall velocity at large distance is ~ 100 
km/s. This hypothetical low density flow is a very rea- 
sonable physical assumption and would explain why most 
spiral galaxies are warped. 
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Accretion due to such a flow is in good accord with the 
observations of the chemical evolution of the Milky Way 
by contributing ~ 1 M Q /yr of low metallicity gas to the 
disc. The High Velocity Clouds, which are presumably the 
accretable material in the Local Group galaxies (Blitz ct 
al. 1999; Braun & Burton 1999; Lopez-Corredoira et al. 
1999; Wakker et al. 1999a) are candidates for a significant 
fraction of the material which fills intergalactic space and 
is accreted. No massive halo is necessary nor high values 
of magnetic fields are necessary, although the presence of 
these elements would not modify qualitatively the present 
conclusions. No calculations in the framework of accretion 
onto the disc are given for a very massive halo but its ef- 
fect would not affect qualitatively the present mechanism, 
as it is explained in §|i]. Models with a very massive halo 
could be generated and different numerical results would 
be obtained depending on the parameters of the halo, al- 
though no major qualitative changes are expected to the 
model presented here since the mechanism of formation 
of warps is dominated by the interaction with the disc; 
only quantitative changes would come from the increasing 
velocity and trajectory curvature of the infalling material 
which would increase the amplitude of the warp, i.e. the 
same amplitude of the warp would be obtained with a 
density even lower than 10~ 25 kg/m 3 . 

Several mechanisms can generate warps but, among 
them, accretion of an intergalactic flow seems to offer a 
very plausible scenario: it is quantitatively consistent with 
many observations and works independently of other in- 
gredients of galaxies and their structure. 
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Appendix A: Numerical calculation of a(R) 

The solution of *^M[ a (JJ)] = 0, from ©, can be 
achieved by means of a numerical calculation with N dis- 
creet values of R, i.e. 

H(a) = 0, (A.l) 
a = (ai, ...,a N ); H = (Hi, H N ), 

ai = a(Ri); Hi = ^(Ri). 

ari 

Newton-Raphson's iterative method is appropriate for 
this kind of numerical calculations. The iteration k + 1 is 
given by 

a k+1 = a k - W- 1 (a k )H(a k ), (A.2) 



Wu ... W 1N 



w N1 ... w NN ) oa i 

By means of this application, we can calculate a(R). 
The first iteration is taken with a(R) — \6q\. 



